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1. Introduction 

One method available to prove the Schauder estimates is Neil Tmdinger's method of molli- 
fication ([Tru|). In the case of second order elliptic equations, the method requires little more 
than mollification and the solid mean value inequality for subharmonic functions. The method 
was lated adapted to the parabolic setting by Xu-Jia Wang in HWanL however in that presentation 
Wang uses an auxiliary estimate coming from the fundamental solution of the heat equation and 
the mean value of property of subsolutions of the heat equations is not used. Our goal in this 
article is show how the mean value property of subsolutions of the heat equation can be used in a 
similar fashion as the solid mean value inequality for subharmonic functions in Trudinger's orig- 
inal elliptic treatment, providing a relatively simple derivation of the interior Schauder estimate 
for second order parabolic equations. 

2. Preliminaries 

For an open subset U G W^, we denote the corresponding open parabolic domain Ux{0,T) C 
denote a (backwards) parabohc cyhnder by Qj? = Bji x {t — R"^ ,t) . We often 
notate a point {x, t) € Ut by X. A (second order parabolic) moUifier is a fixed smooth function 
p £ C^{R'^+^) with //jjd+i pdxdt = 1. For r > we define the scaled moUifier 

1 ( X t 

(2.1) pT{x,t) := —r-j^p 



Let U E M'^+^ and u G Lj^ci^)- For < r < d{X, dP), the mollification of u is given by 

1 f f I X — y t — s\ 
Urix, t) := P[ ~pr ) ^) 

and satisifes spt Ur C Ur, where Ur = {X € U : d{X, dU) > r}. 

The parabolic distance between two points X = {x, t) and Y = (y, s) is defined to be 

d{X,Y) := max{|x - y\, \t - s\'^/'^}. 

We use both sup and [ • [o to denote the supremum of a function. The Holder seminorm is defined 
by 

\u{X)-u{Y)\ 

[uUut ■■= sup 

x^y&Ut a[X,y)" 
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and the Holder norms by 

2 

\u\2,1;Ut '■= ^\dxU\0;UT + \dtu\o;UT 
k=0 

a ; Ut 



The set of functions 

{U € C'^'^{Ut) ■■ [u]2,l,a;UT < 0°} 

endowed with the norm \u\2m,i,a; Ut is called a Holder space. Written out in full the norm is 

2 

\u\2,l,a;UT ■= '^\9x'^\0;Ut + \dtu\o-UT + [t?x^]a;i7T + [dtu]a-UT- 
k=0 

We sometimes use 
e.g, 

[5^'"^'u]a;t/T = [dl'^]a;UT + [dMa-Ur^ 

to compress notation a little. 

We begin by establishing some very basic estimates that will be use throughout. 

Proposition 2.1. We have Ur G C^. 

Proposition 2.2. Let u G L]^^{U). The following estimates hold: 

(2.2) krlo;!/^ < \u\o-U.r 

(2.3) \dl&lur{x,t)\o,u. < CT-'-^^\u\o.Ur- 
Proof. To prove (12.21 ), we have 

1 f f fx — y t — 

^Ho;c/. J J p[—^^-^) dyds 

And for Q: 

dl,d{ur{x,t) = jj^ dldlp (^^-y^,*-^^ u{y,s)dyds 



Proposition 2.3. Let u € Cf^^{U). The following estimates hold: 

(2.4) \Ur{x,t) -u{x,t)\Q.Ur < T'^[u]a;Ur 

(2.5) \dldiur{x,t)\^,ur < Ct'^-'-^^IuU-u^. 



□ 
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Proof. For estimate (12.41 ) we have 

1 If fx — y t — s\ 

Ur{x, t) - U{X, t) = P[ ) '^^ ~ 

< oscu^ u 

To prove the second estimate we have 

dxdiurix,t) = J J dldip(^^—^,*-^^u{y,s)dyds 

The mollifier p is has compact support on Ur and so the last term vanishes by the Divergence 
Theorem. Continuing, we have 

dxd{ur{x,t) = jj^ dl,dip(^--^,^—pf^ {u{y,s) -u{x,t))dyds 

< Ct"*^^-' oscq^ u 

n 

3. The interior elliptic Schauder estimate for functions of compact support 

To motivate things a little for the parabolic setting, we first briefly show how Trudinger's 
method works in the elliptic realm by treating the Poisson equation. The crucial ingredient in 
Trudinger's method is the following norm equivalence: 

Lemma 3.1. Let U C M"^ be an open bounded domain and u G C"{U), where a G (0, 1). Let 
R > and tq be small constants, both of which will be fixed in the course of the proof. There 
exists a constant C = C{d, a) such that the norm equivalence 

-p;[u]a;Bii ^ sup T^"" |(9^-U^ [q; < C[u]a;B^. 

(-^ 0<T<TO 

is valid. 

Proof. The inequality on the right follows directly from equation (12.51 ) (the elliptic version) by 
choosing the appropriate values for the indices i: choosing i = 1 (there is no j in the elliptic 
mollifier) gives 

\dxUr\0;BR < CT"'^[u]a- Br- 

The first inequality requires a little more work. Let x,y G U satisfy |x — ?/| < dx/2 and let r 
satisfy < t < tq < dx/2, where tq will be fixed shortly. For \x — y\ < R = by the 

triangle inequality 

\u{x) - u{y)\ < \u{x) - Ur{x)\ + \ur{x) - Ur{y)\ + \ur{y) - u{y)\ 

< 2t°[u]o;Bb + \dxUr\o;BR\x - y\. 
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Set T = e\x — y\, where e < 1/2. Factoring out and dividing by \x — y[" we find 

Hx) - u{y)\ <\x- yl" (2e"N«;B^ + e°-V^— |a,n,|o;B^) . 

Ciioose e = 4~^/" then set tq = 4~^/". (i2:/2. Taking the supremum over r G (0, tq) completes 
the proof. □ 

We now derive the interior Schauder estimate for Poisson's equation. For simplicity, we only 
consider solutions with compact support in a ball Br. The method extends in the usual way to 
more general equations and domains by using cutoff functions and Simon's absorption lemma 
(see, for example, [iSimJ ). Suppose that u G Cc''^{Bpi) solves 

-a}^ {x)diju{x) = f{x), 

where we assume a*-',/ G C°'{Bji), and there are constants A > and A < oo such that 
< a^-^iiij < We proceed by the method of freezing coefficients, and accordingly 

fix a point xq G M'^ a rewrite the above equation equation as 

-a'\xo)diju{x) - f{xo) = (a'^ix) - a'^ {xo))diju{x) + /(x) - /(xq) 
(3.1) ■.= gix) 

After a coordinate transformation, we can assume without loss of generality that a*-' (xo) = 6^^ , 
that is we can assume (13.11 ) is the Poisson equation. Mollify equation (13.11 ) to get 

-AUr - /(Xo) = gr 

and then differentiate thrice with respect to x to obtain 

-Ad^Ur = dlgr- 

Using inequality (12.51 ) we can estimate 

\dlaTW,BR < C{d)T'^\g\Q.BR 

Recall the solid mean value inequality for subharmonic functions: If v solves —Av{x) < on a 
ball Bji{x) C M.'^, then v satisfies 

. N C(d) f , , , 
v{x) < — / v{y)dy. 

To apply this inequality to our situation, noting A|xp = 2d, we have 

Thus the function d^Ur + \d'^gr\o-^BR\x\^ / {2d) is subharmonic and applying the mean value 
inequality and estimating we obtain 

1 



\dlur{xo)\<C{d), 



dyUriy)dy 



+ R^\dl9rW,Bn 
^OSCBndlurix) + T^^R^'^°' {[a]a-BR\dlu\o. Br + [/IqiBjj)^ 



- ^^^^ ' R 

< C{d)(R''-\dlu]^,BR+r-''R^^''{[aUBR\dluW,BR + 
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We set R = Nt and return to the original coordinates to find 

T'-^\dlur{x^)\ < C(d,A,A) {N''-'[dluU,Bn+N^^''{[a]a;Bn\dlu\o;B^ + , 

then taking the supremum over r and using the norm equivalence we obtain 

[dluUBn < C{d,X,A,a) {N"~^[dluUB^ + N^-^^iiaU-BM^kEn + [f]a;Bn)) ■ 

Choosing sufficiently large and using the Holder space interpolation inequality on the right 
gives the desired estimate, namely 

[dlu]o,-Bii < C{[f]a;Ba + \u\0;Br), 

where C depends on d, A, A, and a. 

4. The interior parabolic Schauder estimate for functions of compact 

SUPPORT 

Having given a feel for Trudinger's method, we move on to use this method to derive the 
Schauder estimates for second order parabolic equations. The crucial equivalence of norms 
lemma in the parabolic setting is the following: 

Lemma 4.1. Let U C M'^+^ be an open bounded domain and u G C'^{U), where a G (0, 1), 
and R > and tq be small constants, both of which will be fixed in the course of the proof. 
There exists a constant C = C{d, a) such that the norm equivalence 

^[u]a;QR< sup { T^"" j^^j^U^ |o; + T^"" l^tU^ |o;Q« } < C[n]c,; . 
C' 0<r<To 

is valid. 

Proof. The second inequality follows directly from equation <\2.5\ by choosing the appropriate 
values for the indices i and j. To prove the spatial part of the second inequality, choosing i = 1 
and j = in estimate (12.51 ) gives 

\dxUr{x,t)\o.Q^ < Ct°"^ Q^. 

The temporal estimate follows similarly. Moving on to the proof of the first inequality, let 
X,Y £ Ut be points such that d{X, Y) < dx /2. To simplify notation a little, set R = dx /2; 
thus d{X, Y) G Qr C Ut. Let r satisfy < t < tq < R, where tq will be fixed shortly. For 

d{X, Y) < R, 

\u{X)-uiY)\ < \u{X)-Ur{X)\ + \Ur{Y)-uiY)\ + \Ur{x,t) -Ur{y,t)\ + \Ur{y,t)-Ur{y,s)\ 

< 2r"[n]a;Q^ + \x - y\\d^Ur\o-Qa + \t - s||(9tn^|o;Q^. 
Set r = ed{X, Y), where e < 1/2. Factoring out d{X, Y)" we have 
\u{X)-u{Y)\ < d{X,Yr {2e''[uU;Q^+e"-\^-''\d.,UrW,Q^+e"-^T^-''\dtUr\o;Qn) ■ 

Choose e = 4~^/" then set tq = 4~^/" • dx/2. Taking the supremum over r G (0, tq) completes 
the proof. □ 
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We now proceed similarly to Poisson's equation to derive the Schauder estimate for the non- 
homongeneous heat equation. Again for simplicity, we only condisider solutions with compact 
support in a parabolic cylinder Qr as the more general estimates follow from this case using 
cutoff functions and Simon's adsorption lemma. Suppose that u € Cc'°'{Qb) solves 

dtu{x,t) - a'-^x,t)diju{x,t) = f{x,t), 

where we assume a*-',/ € C"((5/j)> and there are constants A > and A < oo such that 
A|^P < aJ'-^iiij < Again we freeze coefficients at a point (xq, to) £ Qr and perform a 

coordinate transformation if necessary to get 

dtu{x,t) - Au{x,t) - /(xo,to) = i{a''^{x,t) - a'^ {xo,to))diju{x,t) + f{x,t) - /(xo,to) 

:= g{x,t), 

and then mollify the equation to obtain 

(4.1) dtUr{x,t) - AUr{x,t) = gr{x,t). 

Given the form of the norm equivalence, the desired Schauder estimate will follow if we can 
establish the estimates (for the spatial component of the Schauder estimate) 

(4.2) \d',ur{x,t)\ < C (^loscQ^a2u + i?2|a3^^|o.Q^^ 

(4.3) \dtdlur{x,t)\ < C (^-^oscQj^dlu + R'^\dtdlgr\o;Qn^ , 
and for the temporal part 

(4.4) \dxdtUr{x,t)\ < C (^-^oscQj^dtu + R'^\d^dtgrW,QR 

(4.5) \dtUr{x,t)\ < C oscQj^dtu + R^ldfgrlo^Qi 



We show how to obtain estimate (14.21 ) as the other three estimates follow in the same way. Recall 
the mean value property for subsolutions of the heat equation: If w is a subsolution to the heat 
equation on M'^+^, that is if v solves — Av < 0, then v satisfies 

for each E{x,t;r) C M*^"*"^, where the heat ball E{x,t;r) is the set given by E{x,t;r) = 
{{y,s) e R'^+^ : Ix-yp < Y^-27rslog[r2/(-47rs)], s G {t - /{ATrs),t)}. The radius of 
the heat ball ^ —2-ks log[r2/ (— 47rs)] is often denoted by Rr{s). We refer the interested reader 
to BEval and fEck] for more information on the mean value property and heat balls. Let us now 
show (14.21 ): Differentiate (14.11) thrice in space. Since I^^^Tt-Io; £;|xp/(2n) is independent of time 
we see 

= dlgr-\d^,grkE<0, 
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and hence the function d^Ur + |5^gr!o;£;(xo,to;r)I^P/(2"^) subsolution of the heat equation. 
From the mean value property of subsolutions we have 



(4.6) diuA.,t) < - JJ^^^^^^^^ [^d^,uAy,s) + m^rlo^) J^dyds. 

By translating coordinates we can assume without loss of generality that (x, t) = (0, 0). Estab- 
lishing the required estimates involves evaluation of an integral of the form 

1 f Rr{sY 



ds, 



where a and /3 are given integers. The constants can be computed explicitly, however we are 
only interested in the scaling behaviour with respect to the radius r (and that the integral is 
finite). We compute 



1 1 r ( -2nslog[rV(-47rs)]) 



a/2 



47r 47V 



OO 



Jo 

where the last lines follows from properties of the Gamma function. We point that the integral 
is independent of a and /3 when a = /3 = 2. Returning to (14.61 ). we have 

(4.7) dlur{x,t) <4r"' JJ^d^Ur{y,s)^-^dyds + ^\dlgrW,E jjj^dyds. 
We estimate the first term on the right by 

4r-" / / ^lur{y,s)^-dyds < Cr'^ j\ ( / dlurdy] ds 



r 












r 




47r 


S2 



< Cr-" / ^ ( / oscd^udy ] ds 



C{d) 



r 



The second term on the right of (14.81 ) can be estimated more simply to give 

MdbrkET-'' llj-^dyds < C{dy\d^Mo;E. 

The other three estimates (I4.3l l - (14.51 ) follow in a similar way. For example, to derive estimate 
(I4.3I ). we differentiate (14.11 ) twice in space and once in time. Using the mean value property of 
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subsolutions of the heat equation we obtain 

ff lyP 4r~" ff \y\^ 

(4.8) dtd^Ur{x,t) < Ar'"- J J dsdyUr{y, s)^ dyds + ^^\dtdlgr\o;E J J -^dyds. 

The first part of the estimate follows by integrating by parts in time: 

< — — osce o^u, 

and the second term on right can be estimate in the same way as before to give (I4.3l i. 

The derivation of the paraboUc Schauder estimate now continues in the the same way as for 
the Poisson equation, using the estimates (14.21) - (14.51) . the equivalence of norms lemma and the 
Holder space interpolation inequality. After completing these calculations we obtain obtain the 
desired parabolic Schauder estimate: 

(4-9) [d^Mc.;Q^<C{[fU.,Q^ + \u\o,Q^), 

where C depends on d, A, A, and a. 
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